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A quantum critical endpoint related to a metamagnetic transition causes distinct signatures in the 
thermodynamic quantities of a compound. We argue that, irrespective of the microscopic details of 
the considered material, the diverging differential susceptibility combined with the Ising symmetry 
of the endpoint give rise to a number of characteristic metamagnetic phenomena. In the presence 
of a magnetoelastic coupling, one finds a correspondence of susceptibility, magnetostriction and 
compressibility and, as a result, a pronounced crystal softening, a diverging Griineisen parameter, 
a sign change of thermal expansion a{H), and a minimum in the specific heat coefficient •y{H). We 
illustrate these signatures and their relation on the metamagnetic crossover at 8 T in the prototypical 
heavy- fermion system CeRu2Si2. 

PACS numbers: 75.30.Kz; 71.27.+a; 71.10.HF 



I. INTRODUCTION 

Emergent universality close to critical points is one of 
the most fascinating phenomena in physics. Materials 
with utterly distinct microscopic composition might ex- 
hibit similar behavior close to a second order phase tran- 
sition if they only belong to the same universality class. 
Universality is also expected close to a critical endpoint 
that terminates a line of first-order transitions like, for 
example, in the phase diagram of the liquid-gas transi- 
tion. Such an endpoint is characterized by Ising univer- 
sality, though the Ising order parameter is sometimes not 
simply related to measureable quantities. 

Another example of an Ising critical endpoint can be 
found in metamagnetic materials. Metamagnetism is of- 
ten casually described as a superlinear rise of the magne- 
tization M{H) at some finite critical field Hm- With 
decreasing temperature, such a smooth metamagnetic 
crossover might evolve into a sharp first-order jump of 
M{H) identifying an endpoint {Tcp,Hm) in the phase 
diagram plane defined by temperature, T, and magnetic 
field H}-^ In the context of the metamagnetic material 
Sr3Ru207 it was pointed ouli^'^ that an interesting situa- 
tion arises if the endpoint temperature T^p can be tuned 
towards zero. Top — >■ 0, by a certain external control pa- 
rameter resulting in a quantum critical endpoint (QCEP) 
for Tcp = 0. Such an Ising QCEP differs from its clas- 
sical counterpart at finite Tcp because the dynamics of 
the Ising order parameter has to be taken into account 
explicitly in the quantum case.^ This dynamics then gen- 
erate the temperature dependence in its vicinity giving 
rise to quantum critical scaling. For example, the dif- 
ferential susceptibility at Hm diverges by definition with 



decreasing temperature with a characteristic powerlaw, 
X ~ T~^. Despite the fact that a QCEP can be hardly 
realized in any material due to the required fine-tuning 
Tap = 0, it might nevertheless control thermodynamics in 
an extended temperature and field range if it is only close 
in parameter space. This motivates us to look for uni- 
versal signatures of the metamagnetic QCEP in various 
metamagnetic compounds. 

In the present work, we concentrate on the canoni- 
cal heavy-fermion material CeRu2Si2, which crystallizes 
in the tetragonal ThCr2Si2-structure and shows a pro- 
nounced metamagnetic crossover for field parallel to the 
crystallographic c-direction with a steep rise in the mag- 
netization M (H) at noHm ~ 8 T.^ Over the last 20 years, 
its properties close to the critical field H„i have been in- 
tensively investigated by various experimental methods 
promoting it to be one of the best-studied metamag- 
netic metals, cf. the review Ref. [^. It has been noted 
early on that the metamagnetic signatures are mirrored 
by strong anomalies in dilatometry due to magnetoelas- 
tic coupling ii Close to Hm one observes a remarkable 
large Griineisen parameter a sign change of the thermal 
expansion^ and a strong crystal softening. 1- Moreover, 
the specific heat coefficient ^{H) shows a characteristic 
double peak structure close to the critical field Al- 
though the temperature dependence of thermodynamics 
is anomalous, Fermi liquid behavior is recovered at low- 
est temperatures for all magnetic fields. In particular, the 
differential susceptibility at iJ„i first increases with de- 
creasing T but then starts to saturate at a temperature of 
the order of T* = 0.5 'KM^ In the past, the experimen- 
tal results have been often interpreted within a scaling- 
ansatz for the entropy of the form S{H/ Hm{p), T/To{p)), 
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with the magnetic field H, temperature T, pressure de- 
pendent critical field H^ip) and the temperature scale 
To (p) j^ii^""— This phenomenological approach was quite 
successful to account for the observed relations between 
the iJ-dependence of various thermodynamic quantities 
in the low-temperature limit. However, it did not pro- 
vide an explanation for the huge anomalies themselves 
like, for example, the large Griineisen parameter. Since 
then, a number of microscopic theories based on the peri- 
odic Anderson or Hubbard model have been put forward 
for metamagnetism in heavy-fermion materials^^^— that 
qualitatively explained many of the observed features in 
CeRuaSia. 

Nevertheless, in a previous comparison^^ of CeRu2Si2 
with the ruthenate Sr3Ru207 it was already noted that 
both compounds share similar metamagnetic anomalies 
despite their microscopic differences. This encourages us 
to speculate that some of these anomalies are not specific 
to the microscopies of CeRu2Si2 but are, in fact, associ- 
ated with the emergent universality expected close to a 
metamagnetic QCEP. Note that the importance of crit- 
ical magnetic fluctuations in this material has been an- 
ticipated early on^^ In the following, we demonstrate 
that, actually, some of the most striking metamagnetic 
features in CeRu2Si2 can be naturally explained within 
a QCEP scenario. We argue that the emergent Ising 
symmetry combined with the enhanced differential sus- 
ceptibility give rise to the following generic phenomena 
close to a metamagnetic QCEP: (1) a correspondence 
between susceptibility, magnetostriction and elastic con- 
stants, and, as a consequence, (2) a pronounced crystal 
softening, (3) an enhanced Griineisen parameter, (4) a 
sign change of the thermal expansion, and (5) a minimum 
of the specific heat coefficient, 'y{H), with two accompa- 
nying side peaks. Note, however, that the QCEP itself 
is not realized, neither in Sr3Ru207 nor in CeRu2Si2. 
Whereas in the ruthenate it is masked by a thermody- 
namic phase,— the saturation of the differential suscep- 
tibility in CeRu2Si2 and the concomitant onset of Fermi 
liquid behavior indicates that the QCEP is close but still 
off in parameter space, its distance measured by the sat- 
uration temperature T*. For this reason the sharp in- 
crease in the magnetization in CeRu2Si2 is mostly called 
metamagnetic-like transition or metamagnetic crossover 
in the literature.^ It is an open question whether there 
exists any tuning parameter that lowers the saturation 
temperature T* further to bring CeRu2Si2 closer to quan- 
tum criticality. Attempts to tune the system to the 
QCEP or even to a first order phase transition by ap- 
plying pressure^ or Ge-doping on the Si-site?^ the lat- 
ter being equivalent to negative pressure, failed. In both 
cases the thermodynamic signatures at the metamagnetic 
crossover are even broader compared to the pure system 
at ambient pressure. 

We would like to point out that the concept of a meta- 
magnetic QCEP should be distinguished from a quantum 
critical point (QCP) that separates two phases with dif- 
ferent symmetry. The QCEP is the endpoint of a line 



of first-order transitions that is tuned to zero tempera- 
ture, and as such it is an isolated singular point in the 
phase diagram that is surrounded by a single thermo- 
dynamic phase. In contrast, a QCP separates different 
thermodynamics phases, for example, a paramagnet and 
an antiferromagnet. The latter is actually realized in the 
Ge-doped sister compound CeRu2(Sii_xGea;)2 close to 
concentrations x ~ 0.06 — 0.07.^^ 

For the comparison presented here, we performed new 
measurements on a single crystal of CeRu2Si2 grown 
by Czochralski methodi^ Thermal expansion, magne- 
tostriction and heat capacity measurements were car- 
ried out in a ^He/^He dilution refrigerator mounted in- 
side an 18 T superconducting magnet. A capacitive 
dilatometer^ - made of CuBe and with high sensitivity of 
= 10"^" was used to estimate the thermal expansion 

coefficient ac{H,T) — j^^p- and the magnetostriction 
coefficient Xc{H,T) = xi^ln ^ II c. It is well 
known that in CeRu2Si2 the length change of the sample 
AL shows the same temperature and field dependence 
along all crystallographic axes, but ALc is three times 
larger than ALa-^^ For this reason, we scaled our data 
by a factor of | to get the volume expansion a and mag- 
netostriction A. The heat capacity was measured by a sil- 
ver platform using compensated heat-pulse techniquei^ 
In both experiments we tuned the magnetic field very 
close to the metamagnetic critical field Hm in ultrafine 
field steps of just a few mT. We were able to extend the 
temperature range of our measurements down to 60 mK 
in comparison to former publications. Our results repro- 
duce nicely existing data for the same temperatures and 
fields of a, A^ii^^ii and the specific heatiiii^i^ 

The paper is organized as follows. In Section |ll] we 
discuss qualitatively the characteristic thermodynamics 
expected close to a metamagnetic QCEP, in Section IIIII 
we discuss the metamagnetic signatures and their rela- 
tionship in CeRu2Si2, and we end with a summary in 
section IIVI 



II. UNIVERSAL SIGNATURES CLOSE TO A 
METAMAGNETIC QCEP 

Close to a metamagnetic quantum critical endpoint the 
free energy density, J- — + J- a-, can be separated into 
a background part, To, and a critical part, Jcr, deriving 
from the Ising QCEP 

j-„ = J-„(/i,r,r). (1) 

We assume that yields only a small, featureless back- 
ground contribution to thermodynamics that is sub- 
leading compared to the critical one. The critical part 
(HI) depends on temperature, T, on the magnetic scaling 
field h that is conjugate to the Ising order parameter, 
and a parameter, r, that measures the distance to quan- 
tum criticality. A negative value, r < 0, corresponds to 
a finite endpoint temperature T^p > 0, and for r > 
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FIG. 1: Schematic depiction of planes in the H — T-phase 
diagram with a critical endpoint, which terminates a line of 
a metamagnetic first-order transition. The parameter r is a 
measure of the distance to the QCEP in parameter space, and 
three different cases r^, ri < 0, ri < r2 < 0, ra = are shown. 
For r > the endpoint disappears below the T = axis, which 
corresponds to the situation in CeRu2Si2. The dashed local 
coordinate system close to the endpoint aligns itself with the 
physical coordinates as the endpoint temperature approaches 



zero, Te, 



0. 



only a metamagnetic crossover occurs. The QCEP is 
realized exactly for r = 0. Additionally, it is impor- 
tant how the magnetic scaling field, ft,, is related to the 
physical fields. Usually for endpoints, this relation is not 
evident as, e.g., for the liquid-gas transition, and, as a 
consequence, the Ising symmetry is often hidden. How- 
ever, in our case the situation is more fortunate as we 
are dealing with an endpoint at T = 0. In the limit of 
a vanishing endpoint temperature Top — >■ 0, the line of 
first-order metamagnetic transitions will align itself with 
the temperature axis in the {H, T) phase diagram. This 
follows from the Clausius-Clapeyron relation after tak- 
ing into account that the transition at T = is between 
two ground states of same entropy, see Fig.[T] This align- 
ment has the consequence that the magnetic scaling field, 
/i, of the Ising QCEP can be directly identified with the 
distance to the critical magnetic field 



h\ 



H — H„ 



(2) 



The magnetic scaling field, h, thus controls the distance 
to the QCEP directly on the magnetic field axis and, as a 
result, the Ising symmetry becomes explicit in the (H, T) 
phase diagram. At a finite temperature, there will be a 
superlinear T-correction to h even for the QCEP, r = 0; 
we will later see that for CeRu2Si2 this correction is of 
order 0{T'^). 

Within certain modelsj^i^i^^ the critical free energy ([T]) 
can be calculated and the dependences on its parameters 
can be determined. In the present work, we will follow 
a different route and concentrate on the generic prop- 
erties that all these models for a metamagnetic QCEP 
have in common. We will therefore focus on a qualita- 
tive discussion of the universal metamagnetic signatures 



that derive from two basic assumptions: (i) a diverging 
differential susceptibility x at the QCEP and (ii) its Ising 
universality. 

Probably the most fundamental quantity, that char- 
acterizes the metamagnetic behavior, is the differential 
susceptibility x, 



X 



dm' 



(3) 



At a QCEP, r — Q, the susceptibility diverges upon de- 
creasing temperature, T , or decreasing ft. 



XIqcep 



as T, Iftl 



0. 



(4) 



This serves as our definition of quantum critical metam- 
agnetism. We show below that this increase of x close 
to the QCEP is responsible for all of the striking meta- 
magnetic phenomena. The pronounced increase of the 
susceptibility x in CeRu2Si2 at Hm only saturates at a 
temperature T* = 0.5 K. In the QCEP scenario, this sat- 
uration temperature T* is a measure of a finite positive 
parameter r in Eq. ([1]), ^\h\^[j X ~ I/''- 

Apart from x, we will discuss various other second or- 
der derivatives of the free energy: the specific heat coef- 
ficient 7, thermal expansion a, magnetostriction A and 
compressibility k, defined as 



9^ 



dTdp' 



A 



dHdp' 



(5) 



The pressure dependence of the critical part J-"cr enters 
via the smooth pressure dependence of ft., i.e., the criti- 
cal field, Hra = Hm{p). In the language of field theory, 
the pressure dependence of the other parameter r is less 
relevant and will be neglected in the following. It is con- 
venient to define 



dp 



(6) 



For the small pressures applied in dilatometric experi- 
ments f2„i can be approximated to be constant, flm ~ 
const. 

Because of the scaling form in Eq. ([Ij , the critical con- 
tributions to thermodynamic quantities are not indepen- 
dent. It follows from Eq. (IT|) that the critical part of 
the thermal expansion parallels the T-derivative of the 
magnetization 



' dT 



(7) 



where Mcr — ~dJ-cr/dH . In addition, the critical parts 
of the susceptibility, magnetostriction and compressibil- 
ity are expected to be proportional 



1.1 



(8) 
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Such proportionalities have been observed in 
CeRu2Si2ifi.M and also in SraRusOy^ For a QCEP 
the susceptibility diverges by definition (j4]), which in 
turn implies a divergence in the compressibility k. This 
means that the crystal lattice is destabilized by strong 
metamagnetic fluctuations, and, as a consequence, 
the QCEP is likely to be preempted by a structural 
transition j^^'^^ We interpret the enormous crystal 
softening of up to 50% observed^i^ in CeRu2Si2 as 
a precursor of such a structural instability driven by 
metamagnetic fluctuations. 

A strong increase of Xcr(T) with decreasing T, Eq. 
also has implications for the 7f-dependence of the specific 
heat coefficient, 7cr(i?)- Using the higher-order Maxwell 
relation 

it directly follows from a positive curvature of Xcr (T) that 
7cr(-ff) exhibits a characteristic minimum at Hm- This 
implies that the function 7cr(-ff) must first increase with 
increasing distance from Hm- It is clear that this increase 
of 7cr(-ff) cannot continue indefinitely so that the mini- 
mum at the critical field Hm is likely to be framed by two 
side peaks at finite h. Such a characteristic double peak 
structure is observed in CeRu2Si2i^ and, for thermody- 
namic consistency, we predict it to occur in Sr3Ru2 07 as 
well)^ The saturation of x in CeRu2Si2 at a temperature 
T* is, according to Eq. (0), accompanied by a crossover 
from a minimum to a maximum in 7(i?) close to Hm in 
agreement with experimental observations In addition, 
from the side peaks of the specific heat coefficient 7cr (H) 
we can infer with the help of the relation 

d'fci _ d'^Mcr _ 1 dacr , , 

dH ~ dT^ ~ n.n dT ' ^ ' 

a peak in the temperature dependence of the thermal 
expansion acr(T) at the same finite h. 

Finally, we can exploit the Ising symmetry of the 
QCEP that demands that the free energy is an even func- 
tion of the scaling field h 

T,,{h,T,r)=Tcr-{-h,T,r). (11) 

This Ising symmetry implies, in particular, a sign change 
of the thermal expansion, acr(^) = —o:ci{—h), at /i = 0. 
The position of the sign change in a is a good indicator 
for the strength of the sub-leading temperature depen- 
dence of the magnetic scaling field h{T). If this temper- 
ature dependence can be neglected, the sign change of a 
is located exactly at H ^ Hm and the Ising symmetry 
is explicit in the phase diagram. A finite temperature 
dependence will shift the sign change away from Hm- 
Above, we argued that this T-dependence will be su- 
perlinear due to the Clausius-Clapeyron relation. Char- 
acteristic sign changes in a{H) and peaks in a{T), see 
Eq. (jlO[) . have been observed in the metamagnetic ma- 
terials CeRu2Si2,^ Sr3Ru207^^ and Cai.8Sro.2Ru04.^^ It 



has been shown^S that such sign changes of the thermal 
expansion are generic for quantum criticality and reflect 
the accumulation of entropy. 

In the presence of a QCEP, r = 0, it follows from 
general arguments^ that the Griineisen parameter de- 
fined as the ratio of thermal expansion and specific heat, 
Fcr = c^cr/ijciT), Saturates in the low temperature limit 
to a value given by 

r I ^-it^n ^ fi2) 

^crlqcEP ' H — H 

where the prefactor G is a combination of critical 
exponents4S The sign change in F at Hm reflects again 
the one mentioned above in the thermal expansion. Note 
that if the QCEP is only approximately realized, r > 0, 
as in CeRu2Si2 the divergence Eq. (fT2|) . will be cutoff 
sufficiently close to Hm and the sign change of F will 
not go through infinity but through zero instead^ The 
quantum critical enhancement of F, Eq. (I12p . offers a 
convenient explanation for the anomalous large values of 
the Griineisen parameter observed in CeRu2Si2i^'^ 

At some finite distance to the QCEP, the critical free 
energy ([1]) will be an analytic function of h in the limit 
/i — > 0, and the limiting behavior 

-^cr|T>0 ~ hr{T) - \x.^^AT)h^ + 0{h^), (13) 

is expected. For positive r > 0, this expansion should 
also apply at T = 0. Note that the expansion starts 
quadratically in h due to the Ising symmetry (fTTj). The 
next-order correction of order 0{h'^) is positive, so that 
Xmax(T) corresponds to the maximum value of Xcr(/'') at 
fixed temperature T, see Eq. (j4|). The expression (fT3|) im- 
plies some interesting relations between thermodynamic 
quantities. The susceptibility for vanishing scaling field 
h, Xinax(r) = Xcr{h = 0,r), docs uot Only determine the 
curvature of jcriH) close to Hm, see Eq. but also 
defines the limiting form of the thermal expansion 

acrlT>0 « ^mX[n..iT)h + 0{h') (14) 

where Xmax — f^rXmax- The thermal expansion is ex- 
pected to depend linearly on /i ss — Hm with a pref- 
actor given by the derivative of the critical differential 
susceptibility. Note, however, that we neglect in Eq. 
contributions from the superlinear temperature depen- 
dence of the magnetic scaling field h{T). 



III. COMPARISON WITH CeRuzSiz 

In the following, we discuss in more detail the signa- 
tures in thermal expansion, magnetostriction and specific 
heat of CeRu2Si2 close to the metamagnetic field Hm- 

The phase diagram close to Hm in Fig. [2] summarizes 
the positions of maxima and minima in magnetostriction, 
A(T) and X{H), thermal expansion, a(T) and a{H), and 
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FIG. 2: Magnetic field-temperature {H,T) plane of CeRu2Si2 
close to its metamagnetic field fioHm ~ 7.62 H with the posi- 
tions of extrema in magnetostriction, \{H) and A(r), thermal 
expansion, a{H) and a{T), specific heat ^{H), and the posi- 
tions of vanishing thermal expansion. Note that the drift of 
the positions of maxima in A(_ff), away from Hm spoils the 
explicit Ising symmetry. Critical divergencies in thermody- 
namics are cutoff upon entering the pocket enclosed by the 
positions of extrema in a{H) and a{T) close to Hm', its ex- 
tention can be quantified by the temperature and field scale, 
T* = 0.5K and h* = {Ht,~H*^)/2 = 0.07T//io, respectively, 
which are measures of the distance to the QCEP in parameter 
space, see text. 



specific heat, ^{H) and also shows where the thermal ex- 
pansion becomes zero. For T — > 0, the positions of the 
maxima in magnetostriction \{H) approach and, thus, 
identify the critical field noHm ~ 7.62 T. At finite tem- 
peratures, the positions of these maxima and, similarly, 
the positions of vanishing thermal expansion deviate from 
Hm by a distance proportional to T^, shown by the solid 
and dotted line, respectively. We interpret this deviation 
as arising from a T-dependence of the magnetic scaling 
field h{T) that spoils the explicit Ising symmetry in the 
phase diagram. 

magnetostriction - Magnetostriction as a function of 
field, X{H), has been already presented in Ref . Isl and ITgI. 
In Fig. [21 we focus, alternatively, on the temperature de- 
pendence of our magnetostriction data for different fields 
close to Hm- Upon decreasing temperature, A(T) first 
increases. Away from the critical field, A(r) reaches a 
maximum and then decreases again until it saturates at a 
constant value in the limit T — 0. For fields close to Hm, 
this maximum however disappears and A(T) increases 
monotonously with decreasing T. The magnetostriction 
close to the critical field has the largest absolute values 
identified by the envelope, Xma,x{T), of the set of A(T)- 
curves in Fig. [3) In the following arguments, the pro- 
portionality of magnetostriction and differential suscep- 
tibility X will be of importance; this correspondence was 
impressively demonstrated in Ref.flGl Here, we show this 
correspondence again by comparing the envelope of A(T)- 




7(K) 



FIG. 3: Magnetostriction A(r) for different magnetic fields 
close to fioHm ~ 7.62 T; panels (a) and (b) show A(r) for 
fields H < Hm and H > Hm, respectively, in steps of 0.01 T. 
For comparison, T-dependence of maxima Xmax(r) in the dif- 
ferential susceptibility x{H) from Ref. 'tI are shown identifying 
a proportionality factor Qm ~ 1-5 ±0.1 Tkbar"^, see Eq. ((SJ. 



curves with the temperature dependence of the maxima 
in x(^^)-data taken from Ref. (red squares in Fig. [3]). 
The proportionality factor determines flm as defined in 
Eq. ®. We obtain the value Qm = 1.5 ±0.1 Tkbar'^ 
which is slightly smaller than the value of 2.0 Tkbar"^ 
estimated by pressure experiments^^ The magnetostric- 
tion at the critical field saturates below the characteristic 
temperature T* = 0.5 K, which we identify as a crossover 
from critical to non-critical behavior associated with a fi- 
nite distance to the QCEP in parameter space. As we will 
see below, the inflection point of A(r) located at T* will 
also determine crossover signatures in other thermody- 
namic quantities. 

The positions of the characteristic maxima in A(T) are 
shown in the {H, r)-plane of Fig. [21 Interestingly, the 
maxima in A(T) and x(T) imply, according to the rela- 
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FIG. 4: Thermal expansion as a function of temperature 
for various magnetic fields close to fioHm = 7.62 T. The 
curves are approximately mirror-symmetric due to the emer- 
gent Ising symmetry of the QCEP. 



tion 



9Aci 
dT 



' dT 
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an extremum in the field dependence of the thermal ex- 
pansion, aci{H), at the same positions in the phase di- 
agram, which we will confirm below. Note that charac- 
teristic maxima in the temperature dependence of x(^): 
which are, according to Eq. ([5]), equivalent to the max- 
ima in A(T) of Fig. [31 have been also observed in the 
metamagnetic material Sr3Ru2074i 

thermal expansion - Thermal expansion data of 
CeRu2Si2 has been presented in Refs. and [l^ Al- 
ready in zero field, H = 0, a{T) exhibits a peak that 
shifts to lower temperatures and sharpens with increasing 
H. At the critical field H,m thermal expansion changes 
sign, and the negative peak in a{T) broadens and shifts 
to higher temperatures for increasing fields H > H,n. 
This behavior is reminiscent of quantum critical nietam- 
agnetism as discussed in some detail in Ref. 37. Simi- 
lar behavior of a{T) is observed near the metamagnetic 
field of SraRuzO,^ and Cai.8Sro.2Ru04.^ A closer in- 
spection, however, reveals that the behavior very close 
to the critical field differs qualitatively from the theo- 
retical expectations for a QCEP, see Fig. |4l When the 
peak position has reached a temperature of the order 
of T* = 0.5 K, it does not shift further towards lower 
temperatures upon increasing H, but its height instead 
decreases to zero and reemerges with opposite sign for 
fields H > H,n- The thermal expansion curve a{T) then 
almost recovers its shape but with opposite sign. We in- 
terpret this qualitative change as a crossover from critical 
to non-critical behavior associated with the temperature 
scale T* . For temperatures T < T*, the thermal expan- 
sion has temperature dependence a oq T characteristic 
for a Fermi-liquid. 
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FIG. 5: Thermal expansion as a function of field H for vari- 
ous temperatures obtained by extrapolation from the data set 
in Fig. |4l Panel (a) and (b) show temperatures T < T* and 
T >T*, respectively, with T* = 0.5 K and steps AT = 0.05 K 
and AT = 0.25 K, respectively. Their slope close to the 
critical field is given by dHCi\H=Hm ~ drXmax whose T- 
dependence is shown in the inset (c). Its extremal value for 
the slope, min{9TAmax} ~ —0.001 (T K)~^, at T* determines 
the bundling slope of the Q(ff)-curves close to Hm in panel 
(a). 



The dense data set in Fig. S] allows to discuss the 
magnetic field dependence of the thermal expansion at 
a given temperature, a{H). The extrapolated curves are 
shown in Fig. [S^a) and (b) for temperatures T <T* and 
T > T*. respectively. For low temperatures, T < T*, 
oi{H) has a point refiection symmetry located at H ~ Hm 
and a — which is characteristic for an emergent Ising 
symmetry close to a critical endpoint. The absolute 
value of thermal expansion |a(i/)| first increases upon 
approaching the critical field Hm but after reaching an 
extremum it decreases and vanishes at Hm- The posi- 
tions of these extrema in a{H) are shown in Fig. [21 As 
argued above, these positions in the phase diagram coin- 
cide with the positions of corresponding maxima in A(T). 
In the low temperature limit they extrapolate to the fields 
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T(K) 



FIG. 6: The specific heat coefficient 7(-ff) exhibits a mini- 
mum/maximum crossover close to the critical field Hm with 
decreasing temperature, see inset. The main panel shows the 
curvature djj^{H) (black circles) at the extremum close to 
Hm as a function of T, which is compared to c^l^Amax/fim 
(red squares). By virtue of Eq. ®, the minimum/maximum 
crossover of is determined by an inflection point in the 



specific heat - Specific heat data have been reported in 
Ref. [Tl|. The specific heat coeiiicient 7 is enhanced close 
to the critical field by almost a factor of two compared 
to its zero field value. However, one observes a sharp 
single peak close to Hm only at lowest temperatures, 
T < T*. At elevated temperatures, T > T*, a double 
peak structure in 7(-ff) with a minimum at Hm is found. 
In the following, we focus on the behavior of 7 close to 
the critical field (see inset of Fig. |6]) and, in particular, 
analyze the minimum/maximum crossover at T*. The 
fitted curvatures of j{H) at the extremum close to Hm 
are shown in Fig. [S] for each available temperature. The 
result is compared to the numerical second order deriva- 
tive of Aniax(T) or, equivalently, Xmax(r) of Fig. El As 
anticipated from Eq. ^ , the two quantities agree within 
the error bars. In particular, the minimum/maximum 
crossover is identified with the inflection point of Xmax(T) 
at T « r* where d^XmaAT) = 0. In the phase diagram 
Fig. [2J we also show the positions of the side peaks in 
7(i/) (not shown in Fig. [5]) that coincide within the er- 
ror bars with the positions of the maxima in a(T) as 
expected from Eq. pH)) . 



susceptibility, 



max 



: 0, at T* 



0.5 K. 



IV. SUMMARY 



^oiJ* = 7.55 T and fioH^ = 7.68 T; their distance iden- 
tifies a finite field scale h* = (i7* - iJ* )/2 0.07T/X(7^ 
attributed to the vanishing of the extrema in A(T) and 
a{H). This finite field scale Hoh* = 0.07 T indicates that 
the increase of \a{H) \ while approaching Hm always gives 
way to a decrease even at lowest temperatures. Hence, 
it signifies a crossover from critical to non-critical behav- 
ior at T = and is thus a magnetic field analog of the 
temperature scale T*. 

Furthermore, note that the curves a{H) in Fig. ^a.) 
tend to bundle near iJ,„ to a line with constant slope. 
Generally, the slope of a{H) can be identified with the 
derivative of magnetostriction OtX or, equivalently, of the 
susceptibility dxX- The numerical derivative of XamxiT), 
the envelope in Fig. [31 is shown in Fig. ^c). It is its 
minimum value, min{9TAmax} — —0.001 (T K)^^, at T* 
that identifies the bundling slope in Fig.jSl Moreover, the 
two lines in the phase diagram that identify the extrema 
of a{H) are also the boundaries of the region where the 
analytic expansion of the critical free energy in the scal- 
ing field holds. Well within this region the thermal 
expansion depends linearly on H — Hm as expected from 
Eq. (HID. 

For larger temperatures, T > T*, the thermal expan- 
sion a{H) is shown in Fig. [Sfb). Apparently, the field 
range where a{H) is approximately linear increases, see 
Eq. (|14p . its slope decaying with increasing temperature 
in accordance with the behavior of drXmax- In contrast 
to the low-temperature limit, the curves do not intersect 
at H = Hm and a = anymore, which we attribute to 
the incipient T-dependence of the magnetic scaling field 
h that breaks the explicit Ising symmetry. 



We presented a comprehensive discussion of the uni- 
versal thermodynamic signatures, which emerge close to 
a metamagnetic quantum critical endpoint. We argued 
that (i) the diverging differential susceptibility together 
with (ii) the Ising symmetry of the QCEP account for the 
following characteristics of critical metamagnetism: (1) a 
proportionality between susceptibility, magnetostriction 
and compressibility in the presence of a magnetoelastic 
coupling and, as a result, (2) a pronounced crystal soft- 
ening, (3) an enhanced Griineisen parameter, (4) a sign 
change of the thermal expansion at the critical field, and 
(5) a minimum in the specific heat coefficient ^{H). The 
latter minimum directly follows from the Maxwell rela- 
tion d^x = t^|f7 and the positive curvature of the suscep- 
tibility. A minimum in "f{H), in turn, implies two side 
peaks, and we showed that their positions in the {H, T) 
plane coincide with the extrema in the thermal expan- 
sion a(T) as a function of temperature. The Ising sym- 
metry of the endpoint ensures that the sign change of the 
thermal expansion ot{H) due to entropy accumulation^ 
occurs close to the critical field Hm- 

As an example, we discussed the metamagnetic com- 
pound CeRu2Si2, which shows pronounced metamag- 
netic signatures that saturate, however, at very low tem- 
peratures, T* = 0.5 K, and close to the critical field, 
\H-Hm\ <h* ^ 0.07 T/ig^ We argued that outside this 
regime in the phase diagram behavior for a metamagnetic 
QCEP is expected. We presented new high-precision 
data close to the critical field, that allowed us, in partic- 
ular, to analyze the crossover from critical to non-critical 
behavior associated with the temperature and field scale, 
T* and h* , respectively. We demonstrated that the onset 
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of saturation in the differential susceptibility xi'^)\H=H^ 
at T* leads to an inflection point that accounts for the 
minimum-to-maximum crossover in the specific heat co- 
efRcient "/(H) at Hm- Furthermore, we demonstrated 
that the derivative dTxi'^)\H=H^ determines the tem- 
perature dependence of the thermal expansion close to 
the critical field. The inflection point at T* is reflected 
by an extremum in a{T) at the same temperature as the 
critical field is approached. We also demonstrated that 
the magnetostriction A(T) and, equivalently, the differ- 
ential susceptibility x(r) has a maximum as a function of 
temperature only for fields \H — H„i\ > h* allowing us to 
identify the field scale fJ,oh* — 0.07 T. Such characteris- 
tic maxima have been observed also in the metamagnetic 
compound Sr3Ru207r41 We pointed out that these max- 
ima coincide with thermodynamically equivalent extrema 
in the thermal expansion a{H). 

The focus of the present work are the qualitative ther- 
modynamic signatures close to quantum critical meta- 
magnetism, irrespective of the microscopic details of the 
material and the precise model describing the critical dy- 
namics of the order-parameter fluctuations. We checked, 



however, that all signatures discussed here for CeRu2Si2 
are reproduced qualitatively within the QCEP-model in- 
troduced in Ref.|j, which will be the subject of a separate 
publication. 

The universal metamagnetic signatures have been here 
illustrated on the heavy- fermion compound CeRu2Si2. 
We hope that this work will motivate further experimen- 
tal investigations to identify additional materials that are 
close to a metamagnetic QCEP. 
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